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Summary. — By a length scale analysis, we study the equilibrium interactions be- 
tween two like-charge planes confining neutralising counter-ions. At large Coulombic 
couplings, approaching the two charged bodies leads to an unbinding of counter-ions, 
a situation that is amenable to an exact treatment. This phenomenon is the key to 
attractive effective interactions. A particular effort is made for pedagogy, keeping 
equations and formalism to a minimum. 



1. Introduction 

Consider two identical charged macromolecules alone in an electrolyte, and assume 
that Coulombic forces exclusively are at work. By integrating out the microscopic de- 
grees of freedom (the electrolyte), one obtains the effective pair potential between the 
two macromolecules [1]: is it attractive or repulsive? Of course in vacuum (i.e. without 
the electrolyte), the two bodies always repel, but the presence of the electrolyte compli- 
cates the matter, to such an extent that the answer to the above question has been a 
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Fig. 1 - Sketch of the ge- 
ometry considered, and def- 
inition of relevant parame- 
ters. The two charged par- 
aUel plates at distance d are 
neutralised by counter-ions 
of charge —qe, represented 
by the black dots. The sys- 
tem can be viewed as a mini- 
mal model where like-charge 
attraction may be present, in 
some adequate range of dis- 
tance and coupling. 
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rather controversial subject since at least the 1930s [2, 3, 4, 5]. In essence, all early con- 
tributions pertained to mean-field, a treatment for which it can be rigorously shown that 
the effective potential is repulsive [6, 7, 8]. Yet, it has been understood since the 1980s 
that like-charge attraction is nevertheless possible [9, 10, 11]. Since this phenomenon 
is a signature of non mean-field effects, it is not straightforward to build an intuitive 
picture of the underlying mechanisms. It is our goal here to present a simple argument 
where the prerequisite for like-charge attraction become clear, and which provides exact 
results in some limiting sense to be precised below (namely for short distances and large 
couplings). 

The paper is organised as follows. We first lay out the model in section 2 where 
the contact theorem is also reminded. It is the cornerstone of our analysis. We then 
analyse the limiting case where the two macromolecules are distant (section 3), which 
sheds interesting light on the physics at work when the two plates are close, as developed 
in section 4. Our results are fmrther discussed in section 5 and conclusions are finally 
drawn in section 6. 

2. The model 

We shall concentrate on the simplest setting possible, that where the two macro- 
molecules (colloids) are envisioned as two parallel planar surfaces at a distance d, see 
fig.l. Each colloid/plate is assumed to bear a uniform surface charge distribution ere, 
where e is the elementary charge. The solvent is a structure-less medium of dielectric 
permittivity e, the same as the permittivity of the macromolecules themselves (dielec- 
tric image effects are discarded) and the system is in equilibrium at temperature T. In 
between the colloids, mobile counter-ions of charge —qe ensure global electro-neutrality. 
While the microscopic density of counter- ions is strongly modulated (ions repel), the 
coarse-grained density profile p{z), averaged over a plane at a given position z, is a 
smooth function. For the sake of the argument, the counter- ions are assumed point-like: 
usually, a hard core is required to prevent collapse of oppositely charged small ions, but 
we deal here with a system having one type of small ions only, all of the same charge 
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(no salt), and that is therefore free of divergences. From the permittivity e and temper- 
ature T, we define the so-cahed Bjerrum length is = e^/(efcT) as the distance where 
the repulsion between two elementary charges coincides with thermal energy kT {k is 
Boltzmann's constant, and we use CGS units, so that for water at room temperature, 
one has £b — 7 A). Note that global electro- neutrality imposes the following constraint 
on the counter-ion profile: p{z) dz ~ 2a /q. 

The question is: what is the sign of the pressure P, i.e. the force between the two 
plates per unit surface ? To find the answer, it is essential here to invoke a key result, 
known as the contact theorem [12]. It states that P is, quite intuitively, the sum of two 
contributions: 

(1) P = p{0) kT - — a^e^ 

e 

The first one, p{Q)kT, is repulsive and stems from the collisions of counter-ions with 
the plate. The second, 27r cr^e^/e, is the traditional electrostatic pressure, important 
in capacitor physics. It is attractive, and simply means here that the plate located 
at z = and having surface charge ere sees on its right (z > 0) an integrated charge 
ae — qe p{z)dz = —ere from electro-neutrality, that is thus of opposite sign. The 
contact theorem (1) is exact, and remains true when counter- ion excluded volume is 
accounted for. We note that it implies a strong constraint on the contact density p(Q) 
of an isolated macromolecule. Indeed, when d — >■ oo, the pressure P should vanish so 
that eq. (1) implies p(0) = 27t£b(J^- This result shows that the contact density does 
not depend on the valency q of ions. When d is finite, it is no longer possible to know 
p(0) without an explicit and often complicated statistical mechanics treatment, and the 
ensuing pressure is a non trivial quantity. It however only depends on 2 parameters (a 
reduced distance plus a coupling parameter) and we shall see that for strongly correlated 
systems, meaning for instance that a is large, and for small enough d, p{0) becomes 
simple again and yields an interesting equation of state through (1). Before addressing 
that small d situation, it is informative to discuss in more details the case where the two 
plates are at a large distance from each other. 

3. The large distance limit 

When d is large (meaning d ^ a where a is the lateral correlation length which 
we are about to define), the ionic distribution p(z) around a macromolecule is only 
weakly affected by the presence of the other plate. We have furthermore emphasised 
in the introduction that the relevant situation for like-charge attraction is that where 
non mean-field effects are at work, so that the electrostatic coupling should be large 
enough. To quantify such a coupling, one can assume that all ions are collapsed onto the 
plate (which happens at low enough temperature), in which case the mean surface per 
counter-ion is q/cr, which defines a typical distance a oc \fqJo between ions. Comparing 
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their typical Coulombic repulsion q^e^ jiea) to fcT, we obtain the ratio 

Taking the square, we define the coupling constant S ~ ^ixt^i^a [13, 14]. This parameter 
discriminates weakly coupled cases where S ^ 1 and mean-field holds, from strongly 
coupled ones, where S 3> 1 [13, 15, 14]. Large couplings may be obtained in practice by 
increasing a or q (considering multivalent ions such as Cr^"*", spermidine'^"'" having q — 3, 
or spermine^"*" with q = 4.) 

In the subsequent analysis, we consider S 3> 1, where the counter-ions form a quasi- 
2D layer around the (still here isolated) plate: indeed, they arc then in a configuration 
close to the ground state, which is a hexagonal two dimensional crystal, the so-called 
Wigner crystal [16]. In the ground state (T = or equivalently S = oo), it should be 
noted that any counter-ion feels the electric field created by the plate, while the field due 
to other counter-ions vanishes by symmetry: all ions lie in the same plane. Hence, due 
to thermal excitations and provided d is small enough, the ions can move slightly away 
from the plate and explore a region of size /i, where the energy cost for exciting a given 
ion, to leading order in kT, stems entirely from the plate potential. Remembering that 
the bare plate creates for z > an electrostatic potential —27r ere z/e, we can therefore 
estimate the value of /i by writing 

27r , 1 

(3) qe — ae fi = kT =^ fi 



e 2nq Ib<t 

It is instructive to compare this new length scale fx, called the Gouy length, to both a 
and £b [17]: from the condition S ^ 1, we infer that 



(4) fi <^ a <^ q Ib and more precisely 



a q^li 



The corresponding configuration is shown in fig. 2. We now ascertain the consistency of 
neglecting the field due to other counter-ions on a given one: if the contribution due to 
other ions is accounted for, we obtain that an ion moved from its ground state position 
experiences an additional potential in q^z'^/{ea?). For z ^ ji, this yields kTlB n"^ / oc 
kTS-^/"^ < kT. Here as in most of our discussion, the requirement of a large ^ is 
paramount. (^) 



(^) These argument immediately provide the ion density profile around an isolated charged 
plate, showing that as far as the z coordinate is concerned and to leading order in H, the ions 
decouple from each other, thereby forming an ideal gas in an external gravitational-like field. 
The ensuing profile is simply [18, 19, 20, 21] p{z) oc exp(— 2//i). This expression should become 
asymptotically exact when H — >■ oo -as validated by numerical simulations- entailing that the 
profile should obey the contact theorem, valid at all couplings. This is indeed the case: from 
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Fig. 2 " (Left) Side view of the ionic 
atmosphere in the vicinity of a sin- 
gle charged plate (large d limit). The 
large value of the coupling parameter 
H makes that jj, <^ a. (Right) Front 
view of the system. When T = (i.e. 
H = oo and the Gouy length /i = 0), 
^ the counter-ions are exactly located 
on the hexagonal lattice with con- 
stant a. At finite but large H (the sit- 
uation shown), the counter-ions are 
close to these limiting positions. 



Seen from a large distance by a infinitesimal test charge moving along the z axis as 
shown by the dashed line in fig.2-left, the whole structure (plate-|-ions) does not create any 
electric field since it is electro-neutral [more precisely, the electric field is exponentially 
small for z ^ a, and decays like exp(— z/a)]. However, when the test charge is close 
to the plate (z a), and assuming it does not sit on a counter-ion but at maximum 
possible distance from counter-ions (see the dashed line in fig.2-left), the electric field is 
that of the bare plate, 2'!Taz/e. This innocuous remark plays an important role in the 
argument to follow. 

4. From infinite to small inter-plate distances: the unbinding scenario 

In light of the previous discussion, it is only when d becomes smaller than the lateral 
correlation length a that the ionic structure on a plate starts to distort compared to 
the infinite separation case (this expectation is fully corroborated by a detailed analysis 
[22]). The problem then becomes complex, except when d <C a, see fig. 3. In this short 
separation limit, the ions are confined to a quasi-2D geometry and we need to identify 
the position of a given counter-ion by its projection parallel to the plates r|| , and its 
coordinate z. Both coordinates should be considered separately. Projecting the structure 
onto a plate, thereby eliminating z and focussing onto the set of r|| for all ions, we obtain 
a Wigner crystal very reminiscent of that sketched in fig.2-right, with the only difference 
that the surface density of ions is now doubled, and that the lattice constant is hence 
changed according to a — > a/\/2. On the other hand, the behaviour of the z coordinate 
is more subtle. 



normalisation, it follows that p{z) = 27r^scr^ exp(— z//i) so that p(0) = I-kIbcP' . It should not 
be forgotten though that the ideal gas picture only holds for the z degree of freedom, whereas 
the ions strongly repel in the perpendicular direction (parallel to the plate), where they form a 
crystal. 
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Fig. 3 - The two plates in the 
small separation limit. It is 
assumed here that d <^ a, 
which can still be compatible 
with /i <C rf, see eq. (4). 



A given counter-ion in the slab experiences the sum of the potentials created by the 
plates plus a term due to counter-ions interactions. The first plates-ion contribution 
yields a vanishing electric field (constant potential, say 0) in the symmetric plate case. 
As for the second ion-ion interaction and as was the case for the single-plate analysis 
of section 3, it is a small quantity since all counter- ion basically lie in the same plane, 
with an ensuing small electric field created. To be more specific, this inter-ionic energy 
per particle can be estimated as being of order kTis /a^ oc kT id/ of S^/^; it is thus 
small compared to kT provided that d/a <^ S^^/^ or equivalently d/ ^ <^ a}!^. Under 
that proviso, we can neglect the ion-ion interactions: the ions then experience a uniform 
potential, freely explore the available space < z < d, and consequently adopt a uniform 
density profile along z: p[z) = est. From normalisation J pdz ^ '^cr/q and we necessarily 
have p{z) = 2a /{qd). There is consequently an unbinding oi ions which takes place: when 
d is large, the ions arc confined in a narrow region of extension fi around each plate, see 
section 3. When d <ti a, the counter-ions unbind from the plates, along the z direction, 
as a consequence of the (nearly) cancellation of inter-ionic interactions, and can explore 
the full slab width, of extension d. Note that d can be larger than /x, while being small 
compared to a/S-*^/^, see eq. (4). Still, in the transverse direction , the ions are strongly 
correlated and are essentially frozen on the hexagonal Wigner positions. 

Taking advantage of the unbinding phenomenon, the last ingredient of the analysis is 
to recall the contact theorem which allows to cast the inter-plate pressure in the form 

P 2(7 27rcr2e2 /2u 

(5) T7^ = —; 77^ = 2ttIb<J -7- - 1 

kT qd e kT \ d 

Thus, when d > 2/i, the pressure is negative, and the two like-charge plates attract. On 
the other hand, when d < 2fi, the pressure is repulsive and its positive sign stems from the 
penalising entropy for confining ions in too narrow a slab, which becomes overwhelming 
and leads to the divergence P ^ 2akT/{qd) when d — > 0. It can be pointed out that a 
systematic expansion in inverse powers of the (large) coupling parameter S completely 
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confirms our argument [21, 23], and provides the following correction to the equation of 
state (5) 

(6) = 2^tBcy^ - 1 + - 0.6672... + 0{d^/a^)\ 

kl \ a a ) 

Here, the same requirement as above holds for d/a, namely that dja <^ E^^^^. The 
term d/a in the parenthesis in eq. (6) is therefore a small correction compared to the 
others. Expression (6) shows that the ion-ion interactions that were neglected in deriving 
eq. (5) contribute to a small correction, as expected. In other words, the true density 
profile, leading to (6) via the contact theorem, is not exactly fiat as considered, but 
slightly peaked in the vicinity of the plates(^) [21, 23]. In this respect, the equation of 
state (5) is necessarily a lower bound for the pressure. It is all the more accurate as the 
condition d ^ a/S^/"* is met, and since we may have ^ d, P may take values close to 
the maximally attractive bound — 27rcr^e^/e. This is a key point explaining the stability 
of cement pastes [24], that can be viewed to first approximation as made up of parallel 
charged plates with intercalated multivalent ions. 

5. Discussion 

We now examine in more details some questions raised by our treatment, together 
with a possible generalisation. 

5T. Unbinding scenario and ground-state structure. - The argument put forward in 
section 4 for like-charge attraction hinges upon the unbinding of ions from the near 
vicinity of the plates, which leads to a uniform p{z) profile across the slab. We have 
argued that this could occur at large S provided d <€. but it is clear that such a 

flat profile is far from the ground state structure of the bilayer. Indeed, Earnshaw theorem 
states that an equilibrium position cannot be obtained in a system of point charges under 
the action of electrostatic forces alone [25]. This means here that the counter-ions stick 
to the plates when T = 0, and a more refined study reveals that charges adopt the 
pattern shown in fig. 4 [26, 22]. The corresponding p{z) profile is thus a sum of two d 
peaks at z = and z = d, very non-uniform. . . There is of course no incompatibility 
between our uniform smooth profile and the correct singular ground state. When the 
inter-plate distance d and a are fixed (note that both length scales are temperature 
independent, unlike the Gouy length), and T is decreased to in order to realize the 
ground state, S increases and the condition d <C aj'E}!'^ is violated at some point. There 
is then a crossover between the fiat profile and the strongly peaked one that is reached 
in the ground state. In other words, the analysis of section 4 is valid for large S, but 
nevertheless requires that S ^ [a/d)'^ . By working with suitably rescaled quantities, it is 



(^) these two small peaks are precursors of their ground state counterpart, that develop when 
H 3> {a/d)^, see the discussion in section 5'1. 
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Fig. 4 - Front view of tlie ground-state * 
structure for d <C a. The structure 
is a hexagonal crystal of lattice constant 
a/y/2 (shown by the double-arrow), where 
the open and filled symbols correspond to • 
counter-ions in contact with the plate at 
z = or ai z = d respectively. 



nevertheless possible to obtain expression that are strictly speaking correct for S — >■ oo, 
see below. 

5'2. Back to the failure of mean- field. - It is informative to go back to the mean- field 
treatment, the Poisson-Boltzmann (PB) theory [13], which always predicts repulsion 
between the plates [6, 7, 8]. Yet, the contact theorem that has been used here, is obeyed 
by Poisson-Boltzmann profiles. Furthermore, the mean-field profile also becomes flat for 
small enough d. So, what is going on? The answer simply lies in the fact the length scale 
a has no counter-part in the mean-field approach, where the counter-ions are accounted 
for via their profiles, while discarding their discrete nature. As a consequence, it is only 
for d <C /i that the PB profile becomes fiat [27]. Hence, the equation of state (5) also 
holds for PB theory, but only when d <C /i, a regime where P > 0. 



-0.2 I , 1 , 1 , 1 _ I J -0.2 




Fig. 5. - Rescaled and shifted equation of state as a function of rescaled distance: P — 2/d is 
plotted versus d, with P = P/{kT2-K£BO'^) and d = d/fi. The symbols show the Monte Carlo 
data of Ref [28]: the diamonds are for H = 0.5 (relatively weak-coupling), the circles are for 
H — 100 (moderately strong-coupling) and the squares are for H — 10^ (very strong-coupling). 
The Poisson-Boltzmann (mean-field) result is shown by the upper dotted line. The continuous 
curve displays the locus of points where P = 0. Hence, below this curve, the like-charge plates 
attract. 



This leads us to a second remark. It is often stated that PB theory becomes correct 
when S is small enough (see e.g. [29] for a proof in a system having both co- and 
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counter- ions) . Loosely speaking, this is true, but some phenomena are governed in the 
real system (say in a simulation) by the length scale a, and are therefore necessarily 
missed at PB level, no matter how small S is. This is illustrated in fig. 5, which shows 
pressures obtained in Monte Carlo simulations for three different coupling parameters, 
S = 0.5, S = 100 and S = 10^. The curve for S = 0.5 is close to the PB result shown 
by the dotted line provided d is not too small, sec also fig. 6 below. For S = 0.5, we have 
a fi, and we see that the Monte Carlo data depart from the PB prediction for < 1, 
i.e. d < a. In the rescaled units used in the figure, our equation of state (5) reads 

(7) P^l-^- 

d 

Interestingly, for d <C a and S = 0.5 (upper symbols) we see that our expression is 
perfectly obeyed, since P — 2/d — ?> —1. Indeed, the requirement of a flat profile is met 
for d ^ /i. On the other hand, the PB result is P — 2/d ^ [19], as can be seen 

in the figure, which significantly departs from the simulation data. It can be concluded 
that (5) or (7) are not limited to strongly coupled systems, but also provide the limiting 
small d behaviour at arbitrary coupling S. Yet, we emphasise that (5) or (7) only lead 
to attractive behaviour when S is large (attraction requires having fi <^ a); for S = 100 
or S = 10^ and d > 2/x (d > 2) the pressure does indeed exhibit a negative sign: the 
circles and squares lie below the separatrix P = shown by the continuous curve in the 
lower right corner of fig. 5, see also fig. 6. Moreover, while eq. (5) is fairly well obeyed 
by the Monte Carlo data at S = 100, the improved version (6) fares better and is in 
good agreement with the simulation results, as discussed in detail in [23]. When the 
electrostatic coupling is increased even further, the agreement between simulations and 
eq. (5) becomes excellent, see figures 5 and 6. 

To finish with the scaling form (7), we point out that from dimensional analysis, 
the exact rescaled inter-plate pressure can in general be written as a function of two 
arguments only, i.e. Pexact{d, S). This is true for all values of 2. Working at fixed d and 
letting S — > oo to enforce the strong-coupling limit, we obtain d/a cx dS"^/^ which is 
thus negligible compared to S"^/''. This means that 

(8) lim Pexact[d, S) = i-1. 

This is clearly illustrated in fig. 6. Of course, for larger distances than those in the 
figure, the pressure for 5 = 10^ would start to deviate from the prediction (7), since the 
criterion d/a <C 2^^^ would no longer be met (in the present case, d/a ~ S-^/* corresponds 
to d ~ 20 indicating that deviations would start to appear for already d on the order of 
50 or less). Incidentally, fig. 6 also highlights a) that PB holds for large enough distances 
under weak-coupling (diamonds, 5 ~ 0.5) and b) that for large enough 5, the pressure is 
negative provided d > 2. How large should 5 be is not provided by our argument. It has 
been reported that Sc ~ 12 is the threshold value below which no attraction is possible 
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Fig. 6. - Rescaled equation of state P versus d. The symbols are for the Monte Carlo data of 
Ref [28] and have the same meaning as in fig.5. They correspond to the numerical evaluation of 
Pexact- The Poisson-Boltzmann pressure is again shown by the dotted curve and the thick line 
displays the equation of state (7). 



5'3. Large distance behaviour. - At this point, it is interesting to discuss the phe- 
nomenology for large distances. What is the sign of the pressure? Wc have repeatedly 
emphasised that our argument is limited to small enough distances, so that it is not 
informative for the large d physics. Worse than that, it appears that there is to date no 
reliable answer -be it experimental, numerical or analytical- to the question raised, and 
we are confined here to speculation. A plausible scenario is that when the two plates are 
at a large distance from each other, the ionic profile far enough from the plates is dilute, 
so that the local coupling parameter is small in the interstitial region. Mean-field can 
be expected to hold there [18, 30, 31], from which we might infer that the large distance 
pressure should be positive. Hence, for a given coupling parameter S, we may surmise 
a reentrant behaviour, with a repulsive pressure at both small and large distances, and 
an attractive window in between. This expectation may be dimension-dependant, since 
a recent two dimensional study of the same problem has shed doubts on the re-entrance 
of repulsion at large d [32]. It should be outlined though that the log- form Coulomb 
potential takes in two dimension is scale free, so that the coupling parameter in 2D is 
density independent {f3q'^ where /? is the inverse temperature, which unlike S does not 
depend on the macromolecule charge). This is a notable difference with 3D systems. 

5'4. Asymmetric plates generalisation. - It is straightforward to extend our analysis to 
the situations where the two plates bear unequal surface charges. We leave the derivation 
as an exercise, emphasising that the main difference with the treatment put forward in 
section 4 is that the electric field created by the two plates is now non-vanishing, but 
still uniform. It is given by 2i:[ai — cr2)e/e where Oi denotes the surface charge on plate 
i (plate 1 at z = and plate 2 at z = c?) , from which we define the asymmetry parameter 
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C = (72/ (Ti- We can again, for small inter-plate distance neglect the ion-ion interaction. 
The ions thus experience a linear potential, which leads to an exponential profile in 
exp[— z(l — C)/m]: where /i = /ii is the Gouy length of plate 1. Upon proper normalising 
the counter- ion profile (g J pdz — ui +02), and invoking once more the contact theorem 
which reads 

(9) P = p(0) kT -^ale^^ p{d) kT - ^ al e\ 

we obtain [33, 23, 34] 

This expression reduces to (5) when — > 1, as it should. An interesting check for the 
consistency of the approach is to verify that the same pressure is recovered no matter 
which plate is used for evaluating the pressure in eq. (9). Another interesting and 
exact consequence of the contact theorem (9) is that the interactions between a plate at 
arbitrary charge and a neutral one are always repulsive: — Q leads to P = p{d) kT > 0. 

6. Conclusion 

We considered two interacting parallel charged plates, forming a slab where neutralis- 
ing counter-ions are confined. We have presented an argument to unveil the mechanisms 
behind the like-charge attraction that may ensue at large Coulombic coupling parame- 
ter S, following a mechanical route whereas other approaches often rely on an energy 
route [35, 18, 13]. Our argument hinges upon the fact that at small enough inter-plate 
separation, the counter-ions unbind from the vicinity of the plates. From the contact 
theorem, attraction can only set in when the counter-ion density at contact p(0) is below 
its large distance value, and the unbinding phenomenon is thereby the key to like-charge 
attraction. More precisely, a single counter-ion picture prevails, where due to strong 
lateral repulsion (parallel to the plates), the ions move in the perpendicular z direction 
under the potential created by the plates only, as if they did not interact. This viewpoint 
is fruitful and efficient for the planar geometry considered, but it is misleading in two 
respects: first, it may lead to believe that while the single particle image provides the 
exact leading order pressure when S — > oo, the next correction in the S expansion stems 
from two-body interactions, in a virial-like scheme. This is the essence of the virial ap- 
proach of Refs. [19, 28, 17], which has been shown to be incorrect [21, 23]. Second, the 
single particle picture is in general inappropriate, even to obtain the dominant behaviour 
at large S. This appears most clearly when considering the ionic profile for a single plate 
defining two half spaces with different dielectric constants (work in preparation). 

To dominant order in the coupling parameter S, the fact that there is an underlying 
Wigner crystal of ions has no apparent signature on the interactions (this signature should 
be sought in the next to dominant term [21, 23]). However, the presence of the crystal 
is important for allowing the unbinding of ions and in this respect, the description of 
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strongly coupled charged matter is in the realm of low temperature physics, as discussed 
in the review [20] . It is then instructive to remind the value of the coupling parameter 
where crystallisation occurs, which is on the order of Sc — 3.10^ [36]. If strong coupling 
approaches of the Wigner kind, such as that leading to (6) where restricted to this range of 
coupling, they would be of little interest, since it is in practice difficult to exceed S = 100 
(e.g. trivalent countcr-ions are required in water for already large surface charges such 
that cr£^ is of order 1 [17]). However, it has been shown that the predictions derived 
from the Wigner picture apply "down to" much smaller value than 5c, such as 50 or 
sometimes less [23]. This means that the key point is not the crystal in itself, but the 
existence of a strong correlation hole around each counter-ion. 
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